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We study the correlations of time delays in a model of chaotic resonance scatter- 
ing based on the random matrix approach. Analytical formulae which are valid for 
arbitrary number of open channels and arbitrary coupling strength between reso- 
nances and channels are obtained by the supersymmetry method. We demonstrate 
that the time delay correlation function, though being not a Lorentzian, is charac- 
terized, similar to that of the scattering matrix, by the gap between the cloud of 
complex poles of the S'-matrix and the real energy axis. 



1 Introduction 

The duration of particle collisions is an interesting and important aspect of general scat- 
tering theory which is in a sense complementary to the energy representation ordinarily 
used. A collision is characterized in this case by the time delay of particles in the region 
of interaction. Wigner |jl|] was the first who considered the time on which a monochro- 
matic particle with given angular momentum is delayed during its elastic scattering. He 
established the connection of this time delay to the energy derivative of the scattering 
phase shift. The sharper the energy dependence of the phase shift is the longer is the 
time delay. 

Later on Smith extended the time delay concept on many channel problems intro- 
ducing the time delay matrix 
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in the channel space. Here S is the scattering matrix and the summation index c runs over 
all the M open scattering channels. The matrix Q is hermitian; its diagonal element Q'^'^ 
coincides with the mean duration of collision (time delay) in the c-th entrance channel. 



1 INTRODUCTION 



2 



Generally speaking, the delays are different in different channels c. Taking the trace of 
the Smith matrix, one arrives at the simple weighted-mean characteristic 

Q(E) = ^EQ'' = -i^^ln,etS(E) (2) 

of the duration of collisions. Eq. (|^) is just the many-channel version of the well-known 
simple Wigner formula. (Here and below we set h = 1.) 

The time delay turns out to be an especially pertinent concept for the chaotic reso- 
nance scattering encountered in atomic, molecular and nuclear physics ^, as well as 
in the scattering of electromagnetic microwaves |Q in resonance billiard-like cavities. 
The quantity Q{E), being closely connected to the complex energy spectrum of reso- 
nance states, shows in its energy dependence strong fluctuations around a smooth regular 
variation. The two kinds of variation on different energy scales are naturally decomposed 

Q{E) = {QiE))+QfiiE) , (3) 

with an energy or ensemble averaging. By this, the slow energy dependence of time delay 
is revealed whereas the two-point autocorrelation function 

CQ{E,e) = {QdE+'-)Qf,{E-'-)) = {Q{E+'-)Q{E-'-)) - {Q{E+'-)){Q{E-'-)) (4) 

is used to characterize the time delay fluctuations. 

To the best of our knowledge, the first consideration of these fluctuations has been 
made numerically as well as analytically in and |Q in the framework of rather peculiar 
model of resonance elastic quantum scattering on a leaky surface of constant negative cur- 
vature. The noteworthy feature of this model is that the poles of the scattering amplitude 
turn out to correspond to zeros of the famous Riemann (^-function. The real parts of the 



poles are therefore supposed ||10[ to be chaotically distributed similar to the eigenvalues 
of the Gaussian Unitary Ensemble whereas all their imaginary parts (the widths of reso- 
nances) are the same. The latter very specific property partly deprives the model of its 
interest since actual single-channel widths are known to exhibit quite large fluctuations 
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The width fluctuations are suppressed when many channels are open. In this case 
semiclassical approximation can be as a rule expected to be valid. The semiclassical 
analysis of the time delay in terms of closed periodic orbits is given in [0]. It is in 



particular emphasized there that only the tail of the correlation function (^ corresponding 
to the very large values of e can immediately be related to the (short) periodic orbits. 
Quite opposite, the central peak near the point e = is formed as a result of a strong 
interference of many orbits. Therefore, its width describing the long-time asymptotic 
behaviour of the Fourier transform has no direct connection to the classical escape rate 
and has rather to be calculated on the pure quantum ground. This is in line with the 
results of the analysis |jl3[ of distribution of the resonance widths in the three discs 
scattering problem. 



It is now generally acknowledged that the random matrix theory ||T^ represents a 



suitable and reliable foundation for description of local properties of dynamical quantum 
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chaos |TH|. We therefore use below a random matrix model of chaotic scattering to 
calculate the time delay autocorrelation function. We suppose as usual that the number 
of resonances is asymptotically large and use the powerful supersymmetry technique 
jl6| first applied to chaotic scattering problems in [0 . The number M of the (statistically 
equivalent) scattering channels can be small or large or can even scale with the number of 
resonance states. One can treat the latter two cases [M 0^ 03 as a " semiclassical limit" 



in the matrix model. We show here that the time-delay local fluctuations are governed, 
similar to those of the S'-matrix [|T^, by the gap between the real axis and the upper edge 
of the distribution of resonance energies in the complex energy plane. We compare this 
result with that obtained in the framework of the periodic orbit approach. 

In the next section our statistical matrix model is briefly presented. The connections of 
average time delay with the resonance spectrum and S-matrix fluctuations are elucidated 
in sec. 3. After a short description in sec. 4 of the supersymmetry method which we use 
the main analytical results for the time delay correlation function are given and discussed 
in detail in sec. 5. Some numerical results shedding additional light upon properties of the 
time delay correlations are gathered in sec. 6. We close with a brief summary in sec. 7. 



2 The Resonance Matrix Model 



According to the general scattering theory, the evolution of the A^-level unstable system 
formed on intermediate stage of a resonance collision is described [^, ^ ^ by the 
effective Hamiltonian 

n = H -i-fW, W = VV^ . (5) 

The Hamiltonian (|^) acts within the intrinsic A^-dimensional space but acquires, due to 
the elimination of continuum variables, an antihermitian part. The hermitian matrix H is 
the internal Hamiltonian with a discrete spectrum whereas the rectangular N x M matrix 
V consists of the amplitudes of transitions between A^ internal and M channel states. 
These amplitudes are real in T-invariant theory, so that the matrix W, similar to H, is 
real and symmetric. As usual, we neglect the smooth energy dependence of V and W. 
The dimensionless parameter 7 characterizes the strength of the coupling of the internal 
motion to the continuum. 

The poles of the resonance scattering matrix in the complex energy plane are those of 



the Green's function 20, 21, 22 



g{E) = {E-n)-' . (6) 

They coincide with the eigenvalues £n = E^ — |r„ of the effective Hamiltonian Ti with En 
and r„ being the energy and width of n-th resonance state. It what follows, the properties 
of the spectrum of complex energies £n play the crucial role. 

The intrinsic chaoticity of the internal motion of long-lived intermediate system man- 
ifests itself by chaotic fluctuations in resonance scattering and demands a statistical con- 
sideration. Therefore the random matrix approach extending the well-known [^1], |1^ 
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description of chaotic bounded systems has been worked out in It is usually 

assumed that the hermitian part H of the effective Hamiltonian belongs to the Gaussian 
Orthogonal Ensemble (GOE), 

In the limit N ^ oo eigenvalues of H are situated in the interval [— 2A, 2A] with the density 



given by Wigner's semicircle law. Following p2[, we suggest the transition amplitudes 



also to be statistically independent Gaussian variables, 

(0 = 0, {v:v:^) = ^s'^Xm . (8) 

We will use below the ensemble (|^,|) to calculate the average quantities defined in (|^J|). 



3 Time Delay and Resonance Spectrum 

Since we have neglected a smooth energy dependence of the effective Hamiltonian the 
poles Sn in the lower part of the complex energy plane are the only singularities of the 
resonance scattering matrix. Due to the unitarity condition their complex conjugates S* 
serve as ^-matrix's zeros. These two conditions result in the representation 

det S{E) = n . (9) 

Substituting eq.@ in eq.(0), we come to the important connection 

between the time delay and the trace of the Green's function of the intermediate 
unstable system. The time delay is entirely determined by the spectrum of complex 
energies of this system. The collision duration directly reflects the statistical properties of 
resonances. This is in contrast to the scattering amplitudes S^^ which explicitly depend 
also on the transition amplitudes V^. 

The ensemble averaging of eq.(|lO|) gives 



where m < 1 is the ratio M/N and the function 



satisfies the cubic equati 



{Q{E)) = ^Reg{E) 

g{E) = tX^ (tr g{E)) 
1 rwy .E 



(11) 
(12) 

(13) 



)n [|lg 

giE)- 
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The (unique) solution with a positive real part has to be chosen. It can be seen from 



the consideration given in [jT9| that this real part is close to ■^7rp{E) with p{E) being the 



projection on the real energy axis near the scattering energy E of the density of resonance 
levels in the complex energy plane. 

On the other hand, averaging eq.(|l]) directly, we express (Q) in terms of the two-point 
S'-matrix correlation function [jl^ [T^. In the limit of a large number of statistically 
equivalent channels, M ^ 1, scaling with the number of resonances 



Here 19 



iQ) 



Csie) 



.dCsie] 



de 



iTie) 



+ i 



dT{e) 



s=0 



de 



£=0 



e + iT{e) 

with the two smooth functions defined by 



T{e) ^ K{e) T{e) 



Tie) 



A-fg{e/2) 



2 s{e/2) ' [l + 7g(5/2)r 

and we set = for the sake of simplicity. The quantity 

Cs{0) = T(0) = T 



(14) 



(15) 



(16) 



(17) 



coincides with the transmission coefficient T = 1 — | (S) p. With eq.(|T5|) taken into account 
we obtain from ([T^ ) 

where we have designated r(0) as Tq 



de 



T 



As long as the typical values of the quantity T{e) are small as compared to the pa- 
rameter A characterizing the scale of the smooth e-dependence of the function T{e), the 
two factors on the r.h.s. of eq.(|T3p have quite different energy scales. Only the first fast 
varying factor K{e) describes the local fiuctuations whereas the second one corresponds 
to the joint infiuence of all resonances giving rise to the processes with a very short dura- 
tion. The latter came out from eq. (p^Sj) . The average time delay of a non- monochromatic 
spatially small wave packet caused by the formation of a long-lived intermediate state 
[^, ^ is determined just by the factor K{e) |19| 



(r) = -^ 



dK{e) 



de 







(19) 



£ = 



This implies the connection |26, 19 



(20) 
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4 The Super symmetry Method 



Now we calculate the correlation function (^). Taking into account the relation ([lOD , one 
can cast eq.(Q) into the form 

CQiE,e) = ^Re{(tr^(E + |)tr^t(^_£)) _ (tr ^(i? + |))(tr ^t(^_|))| . (21) 
We also define the normalized quantity 

The terms containing two Green's functions with poles at the same side from the real 
energy axis are omitted in (pT]). We will briefly return to this point later. 

In the limit 7 = 0, when the system gets closed, the correlation function ( pTD becomes 



proportional to the GOE density-density correlation which consists of the singular 
term S{7ipe) and Dyson's smooth function —Y2{7ipe). Coupling to the continuum leads 
to appearing of a new energy scale caused by the decay processes. This scale is defined 
[|T9| by the quantity T{e) from eg . ([IBD . One can anticipate a qualitative changing of 
the correlation function to occur on this scale. For larger distances the influence of the 
antihermitian part should fade away and the asymptotics of Cq for e — > cxd is expected 
to coincide with that of the Dyson's function —Y2{'n'pe). 

To perform the ensemble averaging in (|2T|) we use the modification worked out in [jl9l of 



the supersymmetry technique [0. Using the integral representation of Green's function 
as a multivariate Gaussian integral over commuting and anticommuting variables, one 
gains the possibility to accomplish the averaging exactly. With the help of the Fourier 
transformation in the supermatrix space the integration over initial auxiliary supervectors 
is then carried out. Going along this line, one finally arrives at 

F F N'^ 

{t,g[E + -)tTg\E--)) = -_(str(ar/i)strH2))£ • (23) 
Here the shorthand (. . .)c is used to denote the integral 

(. ••)£ = / d[cx] d[a] exp{-NC{cx, &)}{.. .) (24) 
over two 8x8 supermatrices a and a with the measure defined by the Lagrangian IT^ 



\ 1 1 \ 777, 1 

£,{a,a) = - str a^ — -i? str cr— -str (cr(T) + -str ln(a)+— str ^(l-K/ar/) — -e str (ar^) . (25) 

^ Zi Zi 4r 

The diagonal supermatrices appearing above are equal to 

77 = diag(l, 1,-1, -1,1, 1,-1,-1) 



771 = diag(l, 1, 0, 0, -1, -1, 0, 0) r]2 = diag(0, 0, 1, 1, 0, 0, -1, -1) 
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Here we have set the GOE parameter A equal to one. 



The supermatrix a can be decomposed in the following way |]T7 

a = ToaitTo^ (26) 

where Tq is a transformation from a non-compact manifold whereas the matrix is 
diagonalized by transformations from a compact one. This implies a corresponding de- 
composition of the integrals on the r.h.s. of (0) 

(...)£ = J H(^r) d[aR\d[a]exp{-NCR{aR,a)} J dfiexp{-NC^{aR,To)}{. . .) . (27) 

The Berezinian J-'{aR) depends only on the eigenvalues of ctr; d/i is the invariant measure 
of the manifold of non-compact transformations Tq. At last, the Lagrangian (pSj) is splitted 
into two parts, Cr and given by 

CR{a,a) = |strcr|. — sir aR — |str (cr/jcr) + |str In(o-) , 

(28) 

£^(cT^j,ro) = -lest! {aRTo^r]To) + f str \n{l+jaRTo^r]To) . 

Only the second part depends on the non-compact variables. The first one Cr is 
invariant under a transformation by Tq since it is fully absorbed by an appropriate trans- 
formation of a. One can easily verify that the corresponding Berezinian is equal to unity. 

Since the number of resonances N ^ oo, the integrations over ctr and a can be carried 
out in the saddle-point approximation. At the same time, one has to integrate exactly 
over non-compact variables as long as the number of channels M is finite (m = 0). The 
saddle-point approximation becomes valid for the latter integration when the number M 
also tends to infinity (m is finite). We will consider both cases mentioned. To simplify 
formulae we restrict our further consideration to the center of the GOE spectrum E = 0. 



5 Time Delay Correlation Function 

Let us first consider collisions with a fixed number of channels M. The logarithmic term 
in Cn being proportional to the small ratio m does not influence then the saddle-point 
equations in the {aR, (T)-sector. In particular, the term in (|13D containing this ratio has to 
be omitted. The saddle-point equations are trivially solved in this case and at the point 
E = 

^ = -'^^R^ ' crR = V ■ (29) 

With integrations over aR and a being done, the correlation function (^ reduces to the 
integral 

KQ{e) =2Re j d^sii {naij^ii {Ka2) exp|-7rp£: str ai 2"^^^ ln(l + -TQ;i)| (30) 
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over the invariant measure of the non-compact manifold of To-matrices. Here ai^2 are the 

(31) 



4x4 supermatrices defined in [0, the supermatrix k, = diag(l, 1,-1,-1) and 

47 



^ (1 + 7)^ 

is the transmission coefficient (0) calculated in the limit of m = 0. 



The further calculations go along the line described in details in |]T^ and lead to the 
result 



^1 00 00 

JdXoJ dXi J rfA2/i(Ao, Ai, A2)(2Ao +Ai + A2)^ cos{7rp£(2Ao +A1 + A2)} 



X 



'l-TX,? 



M/2 



where 



/^(Aq; Ai, A2 



;i+TAi)(l+TA2) 

(1 — Ao)Ao|Ai — A2 



(32) 



[(1 + Ai)Ai(l + A2)A2]V2(Ao + Ai)2(Ao+A2) 



The dependence of the function Kq on openness of the unstable system is fully con- 
tained in the last factor in (0). If at least one of the quantities M or T is equal to zero 
the threefold integral reduces to the single one [W\ 



1 °° t 1 

Jdtt(^l- -ln(t + l)^ cos(7rprf) + Jdt (^2 " 2 ^) cos(7rp£t) 



(33) 



= 6{npE) - Y2{npe) 
which is just the normalized GOE density-density correlation function. 

Generally speaking, the threefold integral in (|32|) can be investigated for arbitrary 
number of channels M only numerically using the methods developed in (see the next 
section). However, this integral can be simplified if M becomes large enough. Let the 
number M grow still keeping the ratio m = and the product MT = 2TipTv/ (compare 
with (pO|)) fixed. The quantity Vw is just the limiting value of Fq with T and g calculated 



in the limit m = 0. It coincides with the well-known semi classical Weisskopf estimate p5 
of the correlation length of Ericson fiuctuations. Then 



(1-TAo) 



:i+TAi)(l+TA2) 
and one obtains similar to eq.(Bl) 



M/2 



eic^i-TipTw (2Ao +A1 + A2)}, 



(34) 



Ko(e) 



2 
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+ Jdt e(--^r>^*) (2 - ^ In JZy) cos{npet) . (35) 
2 

This is in close analogy with the consideration of the S-matrix correlation function made 



m 



24 



A new convergency factor appeared in the integrals in (^) as compared to (|3^) where 
only the oscillating cosine cuts the integral in the region of asymptotically large t. This 
makes the function Kq finite for all values of e including zero, so that the 5-function 
is now smeared out. The behaviour of KqIs) is quite different in the regions e -C T^r 
and e ^ Tiy. In the first one it is determined by decays and therefore is sensitive 
to the coupling to the continuum. Quite opposite, for large e the behaviour becomes 
universal since the GOE fluctuations described by the Dyson's function Y2 are restored. 
It is perfectly reasonable since an open system cannot be distinguished from a closed one 
during a small time t -C T^. 

The first 7-sensitive domain is widened when the width Tw grows. In the case of small 
pTiY -C 1 (isolated resonances) it is natural to set aside the contribution of asymptotics 
of the integrand presenting (|35|) in the form 



1 



2 00 

J c/t e"^^^«'* ~ 2 ^^^^^^^ ~ ^) cos(7rp£t) + J dte'^P^"^'^ ~ 2 Ft) ^^^^'^P^^^ ' 
2 

The Lorentzian contribution with the width Tw directly traced to the GOE 5-function 
dominates in the domain e < Tw- The sum of the integrals in the second line is negative 
for all values of e and approaches asjTiiptotically the function Y2 from above. We thus 
come to the conclusion that the correlation function vanishes at some intermediate point 
£0 which can be estimated as 




^o^\— (37) 



using the condition 

1 Tw 



71 p (el+Tl,) 



\y2{peo)\ 



The regime of strongly overlapping resonances, pTyy ^ 1, is the most interesting. In 
this case the main contribution in Kq comes from the region of small t. Therefore, the 
second integral in ( P3| ) can be neglected. Dropping then out the small logarithmic term 
in the first integral and extending its upper limit to infinity, we arrive at 

00 2 2 

KQ{e) ^ f dtte^-^^^'^ cos{n pet) = ^^' ^ . (38) 

Corrections to this result are of higher order with respect to the parameter (pTw)"^- The 
function (R3) is not a Lorentzian at all. Decreasing quadratically in a small vicinity of 
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the point e = 0, it deviates subsequently from a Lorentzian, becomes zero at the point 
e = Tw, reaches a negative minimum and approaches at last zero from below. Just the 
correlation function of such a form with Tw substituted by the classical escape rate was 



conjectured in [|T^] as the limiting classical expression following from the periodic orbit 
picture. However, there is no room for the classical escape rate in the matrix models 
considered here. One can see that the found form has in fact quantum grounds. 

One should return to the exact expressions (^,^) if the ratio m is finite. The reso- 
nances strongly overlap in this case. The saddle-point is now found to be 

To = 1 , a = -ia^^ , ctr = g{s/2) 7] , (39) 

where g is the solution chosen in sec. 5 of the cubic equation (0). The sequential 
saddle-point integrations over an, a and then over the non-compact manifold result in the 
expression 

4 

Koie) = -—^ Re ^ (40) 

where the function V[e) defined in (p!6D is just the one appearing when the S'-matrix 



fluctuations are considered |T^. 



The explicit dependence on e gives rise to a sharp variation of the correlation function 
OD in the vicinity of zero if the typical values |r(£:)| <^ 1 (see eq.(^) and the discussion 



below). As long as the ratio m is small, the quantity is small indeed and we can 
neglect its smooth e-dependence for all £ < Fq ~ Fyi/. Eq.(^D|) is equivalent to eq. ( pSj ) 
within this domain. The asymptotic behaviour for large e also does not change since r(£) 
remains restricted for all e. A small difference can appear only for intermediate values of 
e. 

However, for larger values of m the deviation can become noticeable even near the 
point e = 0. In this case the next term in the power expansion 

r(£)^ro + n£ (41) 



with respect to the smooth e-dependence should be taken into account |]T9|. Because of 
the smoothness, the derivative Fq is small. One can see from eq.(0) that this derivative 
is pure imaginary. The form (pSf) is now reproduced again for sufficiently small e , 



with 



4r2 Vl - 

^Q(^) = (£2 + r2)2 ' 



It has been proven in [|19[ that F^, playing the role of the correlation length of the Ericson 
fluctuations, coincides with the gap between the distribution of resonance energies in the 
complex energy plane and the real energy axis. Therefore we come to the conclusion that 
the properties of fluctuations both of the S'-matrix and time delay are described by the 
same quantity, the gap F^, rather than the classical escape rate. 
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Until now we neglected the "one-sided" contribution 

CQ{e) = {Q)'KQ{e) = 



^Re{(trg(|)tr^(-|))-(tr^(|))(tr^(-|)) 



(44) 



to the correlation function @. As long as m = 0, this contribution is of higher order in 
the parameter N~^. However, this is not the case when the ratio M/N is finite. So one has 
to calculate ( ^i]) explicitly. The well-known replica method ||2^ turns out to be sufficient 
for the latter purpose. Dropping here the corresponding rather cumbersome expressions 
we only note that the function Kqle) is entirely expressed in terms of the slowly varying 
g{^) and varies slowly itself. It has got no pronounced resonance behaviour around the 
point e = and constitutes a smooth background for the correlation function. Its value 
at the point e = is approximately equal to 

so that 

^q(0)/Xq(0) 

The ratio is small under the condition 



1 /TrpTo 

2 V 2iV 



npTo < iV or Tq < 1 



(45) 



implying a clear-cut distinction of the local and global scales . Such a scale separation 
is necessary for matrix models to be valid so far as the fluctuations are concerned. 

The obtained form of the e-dependence of the many-channel correlation function Cq is 
close to that found in for the Gutzwiller's model of single- channel chaotic scattering on 
a space of negative curvature. The same values of all resonance widths and the outcoming 
possibility for resonances to overlap are two specific features of the model which are in 
fact in strong disagreement with properties of the resonance spectra represented by matrix 
models. In particular, the single-channel resonances cannot overlap at all in the latter 
models and their widths fluctuate strongly. That is why our result for M = 1 (see 
below) differs noticeably from the correlation function of ref.[^]. The situation changes 
when the number of channels is large. The width fluctuations diminish with the number M 
of channels growing. Since the time delay depends, according to (0), only on properties 
of the complex energies of resonances and not on the number of channels directly, the 
correlation functions become similar in the two quite different cases compared. 

It is worthy to note that the resonances overlapping strongly suppress the time delay 



fluctuations. Indeed, eq.(|36|) gives for isolated resonances 

1 



whereas 



Kq{0) 



npTw 
1 



> 1 



7r2p2r2 



< 1 



w 



when they overlap. The duration of a collision thus becomes a good deflnite quantity in 
the " quasi classical" limit. 
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6 Numerical results 



Excepting a few limiting cases considered above, further analytical study of (|32D is not 
possible and one has to use numerical methods. However, the threefold integral as it 
stands does not suit for numerical computation. A very convenient substitution of the 
integration variables has been proposed in to overcome all difficulties appearing. 



Following this author we reduce the expression (|32D to the Fourier integral 

oo 

Kgie) = JdtF{t)cos{npet) (46) 



with the Fourier transform F{t) given by a double integral of a smooth function quite 
convenient for the numerical work. The asymptotic behaviour of F{t) can be easily found 
explicitly 

+ tltl ■ (^^) 

For a closed system (T = 0) the Fourier transform F{t) tends to unity in the large- 
t asymptotics. This results in the 5-term in the GOE density-density correlation. A 
singularity still survives even for an open system with one or two decay channels. The 
asymptotics (^Tf ) implies square root or logarithmic divergences correspondingly at the 
point e = in these two cases. 

In Fig. 1 the function Kq{x) versus x = pe is plotted for the case of a single open 
channel. The singular behaviour near zero as well as GOE-like asymptotics are shown. 
The dashed line represents the Dyson's function —Y2{'n'x). The calculation was made for 
the value 7 = 1; only some small domain around zero is sensitive to the choice of 7. The 
correlation function Fig.l has little in common with that found in 0. This discrepancy 
is due to the strong fluctuations of single-channel widths in our model in contrast to 
identical widths of all resonances in Gutzwiller's one. 

For M > 2 the quantity Kq{0) is finite and the correlation function approaches, as the 
number of channels grows, the asymptotics given by (|35D . The Fig. 2 demonstrates this 
for the ratio KQ{e)/KQ{0) in the case of overlapping resonances. In asymptotic regime 



(|38D such a ratio is an universal function of the only variable e/Tw- One can see how the 



exact result (^) gets more and more close to this universal behaviour. 

The Lorentzian peak should dominate the ratio KQ{e) / Kq{0) in the domain e/T]y < 
{7ipTw)'2 ^ 1 when resonances are isolated (see (^7|)). Fig. 3 demonstrates this for two 
values of coupling constant 7. 

As it has been mentioned above, the function Kqle) vanishes at some point Eq. The 
position of this point as the function of the number of channels M at several fixed values 
of 7 is shown in Fig. 4 for three different values of the coupling constant 7. It is clearly 
seen that the square root dependence for isolated resonances (see (|37|)) is replaced by the 
linear one for overlapping ones. 
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7 Summary. 



In this paper we have considered the fluctuations of the characteristic time of colhsions in 
the framework of a random matrix model of resonance chaotic scattering. These fluctua- 
tions are entirely due to the fluctuations of the spectrum of complex resonance energies. 
We calculate analytically the time delay correlation function and investigate its proper- 
ties analytically and numerically for different values of the number of channels and the 
strength of the coupling to the continuum. For any values of these parameters this func- 
tion is far from being a Lorentzian. In particular, it vanishes at some point which plays 
the role of the characteristic correlation length of the fluctuations. In the " quasi classical" 
limit of a large number of strongly overlapping resonances this length is given, similar to 
that of the S-matrix fluctuations, by the gap between the upper edge of the distribution 
of complex energies of resonances and the real energy axis. We do not expect that this 
quantity may be connected to the escape rate appearing in the classical theory of chaotic 
scattering. The latter has been conjectured in pSf to be the semiclassical limit for the 
correlation length in chaotic scattering. 
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Figures 

Fig. 1 The time delay correlation function ( |32D versus x = pe for M = 1 and 7 = 1.0. 
The dotted curve is the Dyson's function —Y2{7ix). 

Fig. 2 Overlapped resonances. The normalized function Kq(£)/_K'q(0) versus x = e/Tw 
for three values of M =5, 10 and 20 (dash-dotted, dashed and dotted curves) and 
7 = 1.0 The solid curve is the asymptotic expression (PH]). 

Fig. 3 Isolated resonances. The normalized correlation function (|36|) for pTyy = 0.1 and 
pT]y = 0.01 (dashed and dotted curves), and Lorentzian (solid curve). 

Fig. 4 The zero pso of Kqle) as function of M for three coupling constants: 7 = 0.01 (*) 
(in this case peo has been blown up by a factor 10), 7 = 0.1 (•) and 7 = 1.0 (o). 
Solid and dashed lines are pTw and the dotted curve is the estimate (p7[). 



